Calculations for Magnetism in SQUIDs at Millikelvin Temperatures 
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Here we present details of a calculation that allows us to extract a surface density of unpaired spins 
from flux vs. temperature experiments performed on field-cooled dc Superconducting QUantum 
Interference Devices (dc SQUIDs). 

PACS numbers: 85.25.Dq, 03.65.Yz, 74.40.+k, 74.25.Ha 



Our recent measurements show an unexpected depen- 
dence of SQUID flux on bath temperature in the mil- 
likelvin temperature range [1]. The flux change scales as 
1/T as temperature is lowered. This behavior has been 
observed in both Al and Nb devices, made both with and 
without a wiring dielectric, and prepared in different fa- 
cilities according to different fabrication recipes [2j . Para- 
magnetic impurities in the materials of the SQUID would 
naturally give rise to such a signature, and we interpret 
the 1/T dependence of the flux through the SQUIDs as 
strong evidence for unpaired spins, most likely in the na- 
tive oxides of the superconductors. 

In order to clarify the source of the temperature- 
dependent flux, we have performed a series of field-cool 
experiments in which a magnetic field Bf c is applied to 
a 350 pH square-washer Nb SQUID (with inner dimen- 
sion 200 /im and outer dimension 1 mm) as it is cooled 
through T c ; the field cool freezes magnetic flux vortices 
into the Nb film, with density a v « Bf c /<&o When 
the device is well below T c , the magnetic field is removed, 
and the SQUID is maintained in a flux- locked loop as it 
is cooled to millikelvin temperatures. In Fig. la we plot 
the flux threading the SQUID as a function of temper- 
ature for eight different values of the cooling field. The 
cooling field strongly affects the temperature-dependent 
flux, enhancing or even reversing the polarity of the ob- 
served signal. In Fig. lb we plot the flux change on 
cooling from 500 mK to 100 mK as a function of the 
cooling field; a linear fit to the data yields a slope of 1.3 
$o/mT. Clearly, vortices contribute significantly to the 
measured temperature-induced flux shift. 

The linear dependence of flux drift on vortex density 
suggests the following interpretation of the data. As the 
superconducting films are cooled through T c in an applied 
field, vortices nucleate and relatively large magnetic fields 
are frozen into the films. The observed signal is due to the 
magnetization of unpaired electron spins on the surface 
of the superconductor in the strong fields produced by 
the trapped vortices. For our Nb thin films we expect 
magnetic fields in the vortex to be of the order of 10 mT. 
This field strength yields characteristic temperatures of 
order 10 mK for single electron spins, compatible with 



the low energy scale seen in the experiments. 

Careful analysis of our field-cool data allows us to ex- 
tract the surface density <j s of spins, a key parameter in 
models of 1/ f flux noise from surface magnetism 0, H, 0] • 
We note that there is negligible direct coupling to the 
SQUID from spins polarized out of the plane of the su- 
perconducting films; this is easily understood from reci- 
procity, as surface magnetic fields due to currents in the 
SQUID have vanishing perpendicular component. How- 
ever, the topology of our device is quite different from 
that of a continuous superconducting washer, due to the 
presence of the vortices. Indeed, the thermal polarization 
of unpaired surface spins in the vortex forces a redis- 
tribution of the circulating supercurrents in the vortex, 
due to the requirement to conserve magnetic flux. These 
currents, in turn, couple strongly to the SQUID loop. 
Calculation of the spin density from the data of Fig. 1 
therefore proceeds in two stages: (1) calculation of the 
coupling between a vortex and the SQUID loop, and (2) 
calculation of the coupling between surface spins and a 
vortex. For both parts of this calculation we need to solve 
for the currents in a superconducting washer with radial 
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FIG. 1: (a) Temperature dependence of the flux threading a 
350 pH Nb/AlOx/Nb SQUID, for different fields B fc applied 
as the device was cooled through the superconducting tran- 
sition, (b) Temperature-induced flux change A<3> on cooling 
from 500 mK to 100 mK, as a function of cooling field Bf c . 
A linear fit to the data yields a slope A&/Bf c = 1.3 <I>o/mT. 
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symmetry (in the case of the SQUID this is a simplifying 
assumption); therefore, we first describe the numerical 
solution to this problem. 

This note is organized as follows. In section I, we de- 
scribe calculation of the current distribution in a thin 2 - 
superconducting washer. In section II, we discuss the ^ 
coupling of a vortex to a SQUID. In section III, we de- x 
scribe the coupling of surface spins to a vortex. Finally, n° 
in section IV, we combine the results of sections II and ^ 
III to extract a surface density of spins from the data of 
Fig. 1. 



I. CURRENT DISTRIBUTION IN A 
SUPERCONDUCTING WASHER 

We consider a thin circular washer with inner radius 
Ti and outer radius r Q . The washer is broken into a set 
of concentric loops with width and turn-to-turn spacing 
w. The vector potential A(x) at radius x and produced 
by a current / through a loop at radius r is given by 
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The differential term dO has magnitude sin# in the di- 
rection parallel to the radius vector of x and cos in the 
tangential direction. Upon integration, only the tangen- 
tial direction is non-zero. The flux enclosed in a loop at 
radius x is given by 
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where we have defined 
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and where K(k) and E(k) are the complete elliptic in- 
tegrals of the first and second kind, respectively. The 
functional dependence of A c is plotted in Fig. 1. 

The self and mutual inductances between loops i at 
radius x and j at radius r are given by 



Mij = ir/iox 



A c ([x + w/2]/r) + A c ([x - w/2]/r) 
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where we have approximately averaged over the width of 
the loop by evaluating A c at x ± w/2. The currents in 
loops 1 to N are described by a vector i; the induced flux 
<j) in the loops is then given by 
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FIG. 2: Plot of normalized magnetic potential A c versus nor- 
malized distance x/r. 



In the case of a superconducting washer with vanishingly 
small London penetration depth Az,, the flux quantiza- 
tion condition 
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requires a constant distribution of fluxes through the 
washer. The current distribution is determined by set- 
ting (j> = $0 for all loops and multiplying by the inverse 
of the mutual inductance matrix: 



(9) 



II. COUPLING OF A VORTEX TO THE SQUID 

A single vortex trapped in the SQUID washer induces 
an effective flux offset in the SQUID. The magnitude of 
this effect can be understood by considering a few sim- 
ple examples. First, for a vortex trapped in an infinite 
superconducting sheet, the phase change between points 
above and below the vortex is 7r, corresponding to a $o/2 
flux shift. By symmetry, this phase change is unchanged 
if the vortex is centered in a superconducting strip of fi- 
nite width. A vortex in a SQUID with a narrow washer 
(u — r )<^r will similarly produce a flux offset of 3>o/2 
in the SQUID if the vortex is placed in the center of the 
wire at (n + r Q )/2. For a vortex placed near the inner 
radius r^, the flux offset will approach 1 <£> , whereas a 
vortex placed near the outer radius r Q will give a flux 
offset approaching $0. For vortices trapped uniformly 
throughout the narrow washer, the average flux offset will 
be $o/2. 

To calculate the coupling between a SQUID and a vor- 
tex trapped in the SQUID washer at radius r, we consider 
a radially symmetric geometry that is amenable to the 
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FIG. 3: Coupling f(r) of a vortex to a SQUID: (Blue) r /n 
= 1.02; (Red) r /vi = 5. The average flux offset for the thin 
washer is (/) = 0.5 $o, whereas for the large washer we find 
(/) = 0.14$ . 



calculation technique described above. From symmetry, 
the flux coupling is independent of angular orientation of 
the vortex around the washer. Thus, the flux coupling 
from one vortex is equivalent to n vortices of magnitude 
&o/n placed around the washer with uniform angular 
spacing and identical radius. As n — > oo, the vortices 
are equivalent to a cut in the washer with a net flux 
through the cut of magnitude $0. This problem can now 
be solved using the matrix inversion technique of eq. [9l 
For the solution, we need to calculate the effective flux 
/ that must be applied to the SQUID so that the total 
circulating current is unchanged when a flux quantum is 
applied to the cut. The flux profile in the loops is thus 



<KO = -/, 
cj>{r') = $ - /, 



r < r; 
r' > r. 



(10) 



The offset flux / is solved using eq. [9] with the constraint 
that the total circulating current ^ i is zero. The solu- 
tion of this problem f(r) is shown in Fig. 3 for both a 
narrow and wide washer. 

By solving for f(r) over the range T{ < r < r Q and 
integrating over the area of the SQUID, we compute the 
average flux (/(r^, r )) coupled to the SQUID from a vor- 
tex with a uniform distribution across the washer: 
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Numerical calculations for a circular washer show that 
(/) = 0.5 $0 for a narrow washer, as expected; for our 
device dimensions (inner radius 100 fim and outer radius 
500 /im), we find a coupling factor (/) = 0.14 <E>o- 

What is the physical source of the offset flux? Vortex 
currents, although localized mostly around the core, fall 



off over a long distance. Part of the current must then 
circulate around the center hole of the SQUID washer, 
giving an effective flux offset to the SQUID loop. 



III. COUPLING OF SURFACE SPINS TO A 
VORTEX 

The total magnetic flux through a vortex is quantized, 
and has contributions from magnetic fields generated by 
the currents circulating around the vortex along with any 
externally applied magnetic field. This external magnetic 
field is typically considered to be zero because of screen- 
ing from the superconductor. However, the surface spins 
located within the vortex are strongly coupled to the vor- 
tex, so their magnetic flux cannot be neglected. As these 
spins become polarized, they couple a flux $> V (T) to the 
vortex; the quantization condition requires a redistribu- 
tion of currents in the vortex. The change in average flux 
coupled to the SQUID is thus = — 0.14 <£ v per vortex 
for the thick washer geometry consider previously. 

To calculate $> V (T) we rely on reciprocity: if we know 
the current distribution in the vortex, we can calculate 
the magnetic fields at all points in space. This allows 
us to calculate the coupling to a spin at an arbitrary 
location; the total flux is obtained by integrating over 
the area of the vortex. We first discuss calculation of the 
current distribution in the vortex, then we describe the 
coupling of the vortex to surface spins. 

We model the vortex as a superconducting washer with 
inner radius n = £ and outer radius r Q » £, where 
£ is the coherence length. The current distribution is 
calculated using the techniques described in section I. In 
this case, however, we have a penetration depth Xl that is 
larger than the inner dimensions of the washer, and flux is 
not constant across the superconductor. Instead, single- 
valuedness of the superconducting wavefunction imposes 
the following quantization condition: 



$0 = j> (A + /i A| J) • dr, 



(12) 



where the integral on the right-hand side is London's flux- 
oid. The new term adds additional elements on the di- 
agonal of the inductance matrix: 
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where Sij is the Kronecker ^-function and d is the film 
thickness. Inversion of the matrix Mlou allows deter- 
mination of the current profile in the vortex. A two- 
dimensional solution of the vortex shows that, for the 
materials parameters we are interested in here, the cur- 
rents are uniform along the z direction to within about 
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FIG. 4: (a) Current density and (b) flux distribution in a 
vortex with £ = 30 nm and A = 100 nm. 

10%. Therefore we may simplify the matrix inversion 
part of the calculation by using 1-D discretization and 
replacing 2\\/dw in eq. [J3]by A/w, where A = 2\\/d 
is the thin-film penetration depth. In Fig. 4 we plot the 
current and flux distributions in the vortex for £ = 30 
nm and A = 100 nm. 

The vortex self-inductance L v is easily determined 



from the total current L, 



in the washer: 



L v $o/I v . 



(14) 



For £ = 30 nm and A = 100 nm, we find a vortex self- 
inductance L v = 0.24 pH. 

Once the currents in the vortex are known, it is 
straightforward to calculate the in-plane B r (r, z) and 
out-of-plane B z (r, z) magnetic fields at all points in space 
by summing over the contributions b ryZ (r, z) of the indi- 
vidual loops of current i. For a current loop with radius 
a centered at r = in the z = plane, we have 




where now we have 
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and where again K(k) and E(k) are the complete elliptic 
integrals of the first and second kind, respectively. 



We now consider magnetic coupling between a vortex 
and a surface spin located at radius r. We first think 
of the spin as an infinitesimal current loop with effective 
sensing area A e ff. The flux coupled to the spin by the 
vortex is thus 



= B(r)A eff , 



(17) 



where B(r) = [B z (r) 2 + ^(r) 2 ] 1 / 2 is the magnitude of 
the field at r. Note that the geometrical coupling of the 
spin to the vortex is unity, since the spin is polarized in 
the direction of the local field. The mutual inductance of 
the spin to the vortex is then 
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The flux coupled to the vortex by a spin at r with 
moment m = \±b tanh (jHBB(r)/2kBT) is thus 

d> s _ = VBL v j-B(r) tanh (^^) • (19) 

The total flux coupled to the vortex by spins with surface 
density a s is therefore 



<$> V (T) = fJLBcr s Lv-^- J 27rrB(r) tanh 
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(20) 
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FIG. 5: Distribution of in-plane B r and out-of-plane B z mag- 
netic fields for a vortex with £ = 30 nm and A = 100 nm. Here 
we consider a film with thickness d — 80 nm and evaluate the 
fields on the surface of the film, assuming a uniform current 
density throughout the thickness of the film. 
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where we have introduced the effective spin polarization 
P e ff(T), defined as 

Peff{T) = ^ / 27TrB{r) tanh (^T^) ^ (21) 

From the known current distribution in the vortex we 
can calculate P e ff(T). Note that in eq. [21] we are also 
integrating the current over the thickness of the washer. 
Because we have a uniform current distribution in z, we 
average the fields over the thickness of the washer. 

IV. ANALYSIS OF FIELD-COOL DATA 

As the SQUID is cooled to millikelvin temperatures, a 
change in temperature yields a change in effective spin 
polarization in the vortex, which in turn induces a change 
in flux coupled to each vortex. The total flux change 
<E>(T) at the SQUID is obtained by summing over all vor- 
tices, taking into account the coupling factor 0.14. We 
find 

*(T) = 0.UA SQ a v fi B a s L v P eff (T), (22) 

where Asq is the area of the SQUID washer, and where 
a v is the density of vortices in the SQUID washer. For 
our field-cool experiments, we have a v w Bf c /$o, where 
Bf c is the magnitude of the cooling field. In these ex- 
periments, we see a linear dependence of flux drift A<1> = 
(100 mK)-$ (500 mK) on vortex density. A linear fit to 
the data of Fig. lb yields a slope A$/Bf c = 1.3$o/mT. 
This slope is related to spin density as follows: 

= 0.14-P tJL B * a L v AP eff , (23) 



where AP eff = P eff (100 mK) - P eff (500 mK). For our 
materials parameters, we find AP e ff = 0.037. Using Asq 
0.96 mm 2 and the slope A§/Bf c = 1.3 $ /mT from 
Fig. lb, we extract a spin density cr s = 5.0 x 10 17 m -2 . 
This density of surface spins is compatible with densities 
considered in recent theoretical models of 1/ f flux noise 
from surface spins @, 0] . 
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